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We show that the optomechanical coupling between an optical cavity mode and the two movable 
cavity mirrors is able to entangle two different macroscopic oscillation modes of the mirrors. This 
continuous variable entanglement is maintained by the light bouncing between the mirrors and 
is robust against thermal noise. In fact, it could be experimentally demonstrated using present 
technology. 
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I. INTRODUCTION 

Entanglement is the most characteristic trait of quantum mechanics [1]. An entangled state of a system consisting 
of two subsystems cannot be described as a product (or a statistical mixtures of products) of the quantum states 
of the two subsystems. In such a state, the system is inseparable and each component does not have properties 
independent of the other components. The nonlocal character of entangled states is at the basis of many paradoxes 
[2], and of the deep difference between the quantum and the classical world. The fundamental role of entanglement 
has been reemphasized in recent years after the discovery that it represents an unvaluable resource for quantum 
information processing [3]. In fact, entanglement is at the basis of secure quantum key distribution schemes [4], of 
quantum teleportation [5], and of the speed-up provided by some quantum algorithms [6]. It is generally believed 
that entanglement can be found only in situations involving a small number of microscopic particles. For example, 
a given amount of entanglement is present between two different spins in the thermal equilibrium state of a system 
of many interacting spins (the so-called thermal or natural entanglement [7]). However, for quantum information 
processing, it is the deterministic generation and manipulation of entanglement which is of paramount importance, 
and in these last years a number of impressive experiments has demonstrated the controlled generation of entangled 
states of two [8], three [9] and four [10] particles. Moreover, since entanglement is one of the distinguishing features 
of the quantum world, it is also fundamental to understand how far it can be extended into the macroscopic domain. 
This is important not only to better establish how the macroscopic classical world emerges from the microscopic one 
ruled by quantum mechanics [11], but also for application purposes. For example, entangled spin-squeezed states 
of atomic samples are known to improve the precision of frequency measurements [12], and the accuracy improves 
with increasing number of entangled atoms. A related question is to establish if and how two macroscopic degrees of 
freedom of two different objects can be entangled. With this respect, a striking achievement has been recently shown 
in [13], where the entanglement between the spin states of two separated Cs gas samples containing about 10^^ atoms 
has been demonstrated. At the same time we proposed a feasible experiment [14] in which even a more macroscopic 
entanglement between the oscillating modes of two mirrors with an effective mass of some milligrams can be generated 
by the radiation pressure of the light bouncing between them (see also [15] for a different and extremely idealized 
model for the preparation of motion entangled states of two cavity mirrors). The continuous variable entanglement 
between two mechanical modes could be used to improve the detection of weak classical forces in optomechanical 
devices as atomic force microscopes or gravitational wave detectors [16,17]. 

In this paper we analyze in more detail and further develop the proposal of [14]. In fact, Ref. [14] restricted to 
the case of identical cavity mirrors, i.e., considered, for each mirror, a single oscillation mode with identical effective 
mass, optomechanical coupling, damping rate and, above all, identical resonance frequency. However, [14] showed that 
the entanglement is present only within a small bandwidth around the mechanical resonance, and since in practice 
two mirrors are never exactly identical, it is important to establish the conditions under which entanglement can be 
generated between two mechanical modes with different resonance frequencies, and its dependence on the frequency 
mismatch. 

In Section II we describe the optomechanical system under study in terms of quantum Langevin equations. In 
Section HI we solve the dynamics of the system in the frequency domain, and then we characterize in detail the 
entanglement between the two mirrors. Section IV is for concluding remarks. 
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II. THE SYSTEM 



We consider an optical ring cavity in which two perfectly reflecting mirrors can both oscillate under the effect of 
the radiation pressure force (sec Fig. 1). The motion of each mirror is the result of the excitation of many oscillation 
modes, either external [18,19] or internal [20,21]. The former is important for suspended mirrors since the excitation 
of pendulum modes of the suspension system leads to global displacements of the mirror. The latter corresponds to 
deformations of the mirror surface due to the excitation of internal acoustic modes of the substrate. These various 
degrees of freedom have however different resonance frequencies and one can select the mechanical response of a single 
particular mode by using a bandpass filter in the detection circuit [22]. For this reason we shall consider a single 
mechanical mock^ for each mirror, which will be therefore described as a simple harmonic oscillator. Since we shall 
consider two mirrors with similar design, the two modes will be characterized by different, but quite close, values for 
the frequencies, Oi and f22, and for the effective masses, mi and m2- 

The optomechanical coupling between the mirrors and the cavity field is realized by the radiation pressure. The 
electromagnetic field exerts a force on a movable mirror which is proportional to the intensity of the field, which, 
at the same time, is phase-shifted by a quantity proportional to the the mirror displacement from the equilibrium 
position. In the adiabatic limit in which the mirror frequency is miich smaller than the cavity free spectral range 
c/(2-\/2L) (L is the diagonal of the square optical path in the cavity, see Fig. 1), one can focus on one cavity mode 
only (with annihilation operator b and frequency Wb), because photon scattering into other modes can be neglected 
[23]. One gets the following Hamiltonian [24] 



=1 



-?i6'f6^(-l)^Gjgj +^^76 (/3"e-''^''»*6l'-/3"*e*'^'"'*6) , 

where g, and Pi are the dimensionless position and momentum operators of the mirrors with [qi,Pj] = iSij, Gj = 
{u)b/2L)y/h/mjflj {j = 1,2) are the optomechanical coupling constants, and the last terms in Eq. 1 describe the 
laser driving the cavity mode, characterized by a frequency uto and a power P^" = hwbo\(3^"'\'^ (76 is the cavity mode 
linewidth) . 

A detailed analysis of the problem, however, must include photon losses, and the thermal noise on the mirrors. It 
means that the interaction of the optical mode with its reservoir and the effect of thermal fluctuations on the two 
mirrors, not considered in Hamiltonian (1), must be added. This can be accomplished in the standard way [25,26]. We 
neglect instead all the technical sources of noise, i.e., we shall assume that the driving laser is stabilized in intensity 
and frequency, also because recent experiments have shown that classical laser noise can be made negligible in the 
relevant frequency range [19,20]. The full quantum dynamics of the system can be exactly described by the following 
nonlinear Langevin equations (in the interaction picture with respect to hubb^b) 

b = i{ubo - uJb)b - ib{Gm - G2q2) -^b + ^ (6™ + A„) , 

qj = n^pj , (2) 
Pj = -^jQj + i-yGjb^b - TjPj + , 

where Tj (j = 1,2) are the mechanical damping rates of the mechanical modes, b^"{t) represent the vacuum white 
noise operator at the cavity input [25], and the Langevin noise operators for the quantum Brownian motion of the 
mirrors are £,j{t). The non- vanishing noise correlations are 

{b''^{t)b''^\t')) =S{t-t'), (3) 

i^jmkit')} = Sj,k ckv^ ^coth l^-^j COS [iv{t - t')] - i sin [uj{t - t') 

where fcs is the Boltzmann constant and T the equilibrium temperature (the two mirrors are considered in equilibrium 
with their respective bath at the same temperature). Notice that the used approach for the Brownian motion is 
quantum mechanical consistent at every temperature [26]. 

We consider the situation when the driving field is very intense. Under this condition, the system is characterized 
by a semiclassical steady state with the internal cavity mode in a coherent state and a displaced equilibrium 
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position for the mirrors. The steady state values are obtained by taking the expectation values of Eqs. (2), factorizing 
them and setting all the time derivatives to zero. One gets 

(g,)., = (-VG,-|/3|Vf^,-, 

{Pj)ss = 0, (4) 
/3=(6)«, = V7^/3™/(76/2-iA6), 

where = Ubo — — Gi{qi)ss + G2{qi)ss, is the cavity mode detuning. 

Under these semiclassical conditions, the dynamics is well described by linearizing Eqs. (2) aroimd the steady state. 
If we now use the same symbols for the operators describing the quantum fluctuations around the steady state, we 
get the following linearized quantum Langevin equations 

h = lAbb - iPiGiqi - Gags) -^b + ^b'^ , 

Qj = Qjpj , (5) 
p. = _n .g. + {-yGj{f3*b + /36t) - Tjpj + . 



III. ENTANGLEMENT CHARACTERIZATION 



The time evolution of the system can be easily obtained by solving the linear quantum Langevin equations (5). 
However, as it happens in quantum optics for squeezing (see for example [25]), it is more convenient to study the 
system dynamics in the frequency domain. In fact, it is possible that, due to the effect of damping, and thermal and 
quantum noise, the two mechanical modes of the mirrors are never entangled in time, i.e., there is no time instant in 
which the reduced state of the two mechanical modes is entangled, unless appropriate (but difficult to prepare) initial 
conditions of the whole system are considered. Entanglement can be instead always present at a given frequency. 
In fact, the two mirrors constitute, for each frequency, a continuous variable bipartite system which, in a given 
frequency bandwidth, can be in an entangled state. The Fourier analysis refers to the quantum fluctuations around 
the semiclassical steady state discussed in the preceding Section, and the eventual entanglement found at a given 
frequency would refer to a stationary state of the corresponding spectral modes, maintained by the radiation mode, 
and which decays only when the radiation is turned off. The spectral analysis is more convenient also because in such 
systems the dynamics is experimentally better studied in frequency rather than in time. The same kind of spectral 
analysis of the nonlocal properties of a bipartite continuous variable system has been already applied in Ref. [27] which 
demonstrated the EPR nonlocality between two optical beams of a nondegenerate parametric amplifier, following the 
suggestion of [28]. 

Performing the Fourier transform of Eqs. (5), one easily gets for the mechanical modes operators (j = 1, 2) 



qj{u) = Bj{w)bin{w) + B*{-u)bl{-u) + %i(a;)a('^) + S,- 2(a;)6(a;) 



(6) 
(7) 



where 



^j,fe(w) 



1 



I 



^3-jX3-j(w) 
1 



-iGs-jGs-kWl' 
1 



(8a) 



(8b) 



Olll2Xl('^)X2(w) 



+ 



G 



1 



-i{Ab+ Lv) 



+ i{Ab-Lu)_ 



and Xj{^) = [^ 



UJ 



V*{u) = V{-uj), but B*{uj) ^ Bi-uj). 



(8c) 

iujTj]~^ is the mechanical susceptibility of mode j. Notice that ^.(a;) = Sj^;j(— w) and 



3 



The simplest way to establish the parameter region where the oscillation modes of the two cavity mirrors are 
entangled is to use one of the sufficient criteria for entanglement of continuous variable systems already existing in 
the literature. These criteria are inequalities which have to be satisfied by the product [14,29,30] or the sum [31,32] 
of variances of appropriate linear combinations of the rcscalcd position and momentum operators of the continuous 
variable systems. These criteria are usually formulated in terms of Heisenberg operators at the same time instant, 
satisfying the usual commutation relations [qj{t),Pk{t)] = iSjk [29,31,32], but they can be adapted to their Fourier 
transform, as long as the commutators between the frequency-dependent continuous variable operators are still a 
c- number [14]. This condition is satisfied in the present case thanks to the linearity of the Fourier transform and to 
the linear dynamics of the fluctuations (see Eqs. (5)), implying that the commutators are always c-number frequency- 
dependent functions. 

The paradigmatic entangled state for continuous variable systems is the state considered by Einstein, Podolski and 
Rosen in their famous paper [2], i.e., the simultaneous eigenstate of the relative distance qi — and of the total 
momentum pi + p2. In an entangled state of this kind, the variances of these two operators are both small and it is 
therefore natural to use them. Defining u = qi — q2 and = pi + p2 , an inseparability criterion for the sum of the 
variances in the case of arbitrary c-number commutators is [31] 

{Auf) + {{Avf)<2\{[q„p,]}\\ (9) 

while that for the product of variances is [14,29,30] 

{iAuf)(iAvf)<\{[q,,p,])\\ (10) 

It is easy to see that the condition (9) implies condition (10), which means that the product criterion (10) is easier 
to satisfy, and for this reason we shall consider only the latter from now on. Furthermore, the product criterion (10) 
allows us to establish a connection with Refs. [28], which showed that when the inequality 

{Auf)(iAvf)<\\{[quPi])\', (11) 

is satisfied, an EPR-like paradox arises [2,33], based on the inconsistency between quantum mechanics and local 
realism, which has been then experimentally confirmed in [27]. The sufficient condition for inseparability of Eq. (10) 
is weaker than condition (11), but this is not surprising, since entangled states are only a necessary condition for the 
realization of an EPR-like paradox (see however the recent paper [30] where it is shown that the weaker inseparability 
sufficient condition (10) can be considered as a marker of the existence of generalized, weaker, EPR correlations). 

In order to apply the inseparability criterion (10) in the frequency domain, we have to make the frequency dependent 
operators gj(w) a,ndpj{u)) Hermitian, i.e., to consider the Hermitian component 

for any operator 0{lo). Using the fact that {qj{uj)) = {pj{Lu)) = 0, j = 1,2 and Vw because they are associated to 
fluctuations around the semiclassical steady state, Eq. (10) therefore becomes 

(7^^,_,,)(7^^^+p,)<|([7^,,,7^p,])|^ (13) 

which suggests the following definition of degree of entanglement for the mechanical oscillation modes at frequency u 
of the two cavity mirrors [14] 

which is a marker of entanglement whenever £^(0;) < 1. 

Using Eqs. (8) it is possible to derive the analytic expression of E{u), which is however very cumbersome. The two 
variances in the numerator of (14) are 

iK-J = i {1^1 (^) - ^2(^)1' + m-co) - B2{-iv)f 

+AAi(a;)|Si,i(w) - S2,i(u;)|2 + AA2(u;)|Si,2(u;) - ^2,2{^)\''} , (15) 
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+Afi{uj) [^i4(w)^24(-'^) + ^i,i(-w)^2,l('^)] 

+Af2{uj) [Si^2(w)S2,2(-w) + Si,2(-w)S2,2(w)]} , 



(16) 



with Njiuj) = a;(rj/f2j) co\h(Uio j 2k bT), while the commutator in the denominator of (14) is given by 




(17) 



In Figs. 2-4 we have studied the behaviour of £^(0;) as a function of frequency and temperature, for different values 
of the difference between the two resonance frequencies of the mechanical modes, fii — f72. This is an important 
parameter because wc have seen in [14] that in the case of identical mirrors, the two mechanical modes are entangled 
only within a small bandwidth around the mechanical resonance. Since in practice the two mirrors will never be 
exactly identical, it is important to establish if the macroscopic entanglement is able to tolerate a certain amount 
of frequency mismatch. For the other parameter values we have considered an experimental situation comparable to 
that of Refs. [20,22,34], where the studied mirror oscillation mode is a Gaussian acoustic mode. We have therefore 
considered a cavity driven by a laser working at A = 810 nm and power P^"' = 1 W, with length L = \ mm, detuning 
Ab = 6 MHz, optical finesse T = 25000, yielding a cavity decay rate 7;, = 6 MHz. The mechanical modes have been 
taken with (effective mass mi = TO2 = 23 mg, damping rates Fi = r2 = 1 Hz, and ili = 1 MHz, while we have changed 
the values of 5^2 around those of Qi. These choices yield for the optomechanical couplings Gi ~ G2 — 2.5 Hz. 

Fig. 2 shows E{(jj,T) for no frequency mismatch, ili = ^l2, Fig. 3 refers to the case with f22 — f^i = 10 Hz, and 
Fig. 4 refers to the case with ^2 — ^1 — 20 Hz. In all cases, the region of the lu,T plane where the two mechanical 
modes are entangled is centered in the middle of the two mechanical resonances, i.e., E{lo,T) always achieves its 
minimum at w = (fii + f22)/2. The frequency bandwidth of the entanglement region rapidly decreases with increasing 
temperature, so that, with the chosen parameter values, entanglement disappears above T ~ 4 K. As expected, the 
oj,T region where the two mirrors are entangled becomes smaller for increasing frequency mismatch (compare the 
three figures). Nonetheless these results are extremely interesting because they clearly show the possibility to entangle 
two macroscopic oscillators (with an effective mass of 23 mg) in a stationary way, using present technology. In fact, 
the two modes are still clearly entangled at T = 2 K and with ^2 — ^1 = 10 Hz (ten times larger than the width of 
the mechanical resonance peaks, see Fig. 3), while one is forced to go below T = 2 K when the frequency mismatch 
is equal to 20 Hz (see Fig. 4). 

Differently from temperature and frequency mismatch, and as it can be seen from the involved analytical expression 
above, it is difficult to determine how the degree of entanglement depends upon the other parameters. It can only 
be verified that, as expected, entanglement improves with increasing mechanical quality factor Qj = flj/Tj and 
that it strongly improves with increasing the effective optomechanical coupling constant, which is given by (3Gj (see 
Eqs. (5)). This shows that for achieving even a more macroscopic entanglement (i.e., larger masses), one has to use 
smaller cavities and, above all, larger optical power. The fundamental importance of the effective coupling constant 
(3Gj also helps us to show which kind of entangled state of the two mirrors is generated by the radiation pressure. In 
fact, when the cavity mode intensity becomes larger and larger, the optomechanical interaction tends to project the 
two mechanical modes onto an approximate eigenstate of Giqi — G2q2 (see Eqs. (1) and (5)), which, since in our case 
it is Gi ~ G2, is essentially equivalent to the relative distance qi — q2- The two oscillators occupy a state that, like 
a standard EPR state, has a very small variance of the relative distance u = qi — q2- On the other hand, since the 
radiation pressure docs not have analogous effects on the total momentum w = pi +P2, the state of the mirrors does 
not exhibit such a small value for the variance ^(Au)^^ as the standard EPR state does. Nonetheless, at large optical 

intensities, as shown by the product criterion of Eq. (10), the effect of the radiation pressure force on the relative 
distance is sufficient to entangle the two macroscopic oscillator modes. Moreover, as it can be seen from Figs. 2-4, the 
degree of entanglement E{u) lies even below 1/4 at sufficiently low temperatures, allowing therefore in principle also 
an experimental test of EPR nonlocality with macroscopic oscillators, on the basis of the inequality (11) of Refs. [28]. 
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IV. CONCLUSIONS 



We have shown how the optomechanical coupling realized by the radiation pressure of an optical mode of a ring 
cavity is able to entangle two macroscopic collective oscillation modes of two cavity mirrors. Using parameter values 
corresponding to already performed experiments involving an optical cavity mode coupled to an acoustic mode of the 
mirror (with an effective mass of many milligrams) we have shown that an appreciable entanglement is achievable at 
temperatures of some Kelvin. This continuous variable entanglement is established at a given frequency, between the 
spectrally decomposed oscillation modes of the two mirrors (see also Refs. [27,28] for an analogous spectral analysis 
of the nonlocal properties of the beams of a nondegenerate optical amplifier). One has a stationary entanglement, 
which is maintained by the strongly driven cavity mode as long as it is turned on. Using the degree of entanglement 
E{u)) of Eq. (14), suggested by the inseparability condition of Eq. (10), we have seen that the entanglement is more 
robust when the two mechanical resonance frequencies are equal (Fig. 2), but that it tolerates a resonance frequency 
mismatch of tens of Hz, miich larger than the width of the resonance peaks. The best entanglement is always achieved 
in the middle of the two mechanical resonances and the frequency bandwidth of the entanglement parameter region 
rapidly decreases with decreasing optomechanical coupling and increasing temperatures. 

This continuous variable entanglement between two macroscopic collective degrees of freedom can be experimentally 
measured using for example the three-cavity scheme described in detail in [14]. In such a scheme, a ring cavity is 
supplemented with two other external cavities, each measuring the spectral components qj{ui) and Pj{u)) of each 
mirror oscillation mode via homodyne detection. With these measurements, it is possible to obtain both variances 
i'R-q^-q^) and {TZp^_^_p^). As it has been verified in [14], if the driving power of the meter cavities is much smaller 
than the driving power of the "entangler" cavity mode, the two additional cavities do not significantly modify the 
entanglement dynamics. A simplified detection scheme, involving less than three cavities is currently investigated. In 
fact, the homodyne detection of the entangler mode b provides direct information on the relative distance between the 
mirror modes q\—q2- The measurement of the total momentum quadrature p\ +p2 could be then achieved using the 
result of this homodyne detection and that of the homodyne measurement of the motion of a single mirror provided 
by a second "meter" cavity mode. It is however possible that an even simpler detection scheme exists, using the 
entangler cavity mode only. 

Another important aspect which has to be taken into account is that the motion of each mirror is the superposition 
of many oscillation modes with different resonance frequencies. We can safely verify the entanglement between the two 
considered oscillation modes provided that the other modes of the two mirrors are sufficiently far away in frequency 
so that their contribution at the analysed frequencies is negligible. Moreover, the above analysis also applies, almost 
unmodified, to the case when the two modes belong to the same mirror. 

The possibility to prepare entangled state of two macroscopic degrees of freedom is not only conceptually important 
for better understanding the relation between the macroscopic world ruled by classical mechanics and the quantum 
mechanical microscopic substrate, but it may also prove to be useful for some applications. For example, Ref. [16] has 
showed that entangled states of the kind studied here could improve the detection of weak mechanical forces acting 
on the mirrors as those due to gravitational waves [35]. 



[1] E. Schrodinger, Proc. Cambridge Philos. Soc. 31, 555 (1935). 

[2] A. Einstein, B. Podolski, and N. Rosen, Phys. Rev. 47, 777 (1935). 

[3] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, (Cambridge University Press, Cam- 
bridge, 2000). 
[4] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

[5] C. H. Bennett et al, Phys. Rev. Lett. 70, 1895 (1993); D. Bouwmeester et al, Nature (London) 390, 575 (1997); D. Boschi 
et al, Phys. Rev. Lett. 80, 1121 (1998); A. Furusawa et al, Science, 282 706 (1998); D. Vitali et al, Phys. Rev. Lett 84, 
445 (2000); T. C. Zhang et al, e-print quant-ph/0207076; W. B. Bowen et al, e-print quant-ph/0207179. 

[6] A. K. Ekert and R. Josza, Rev. Mod. Phys. 68, 3733, (1997). 

[7] X. Wang, and P. Zanardi, e-print quant-ph/0202108 and references therein. 

[8] E. Hagley et al., Phys. Rev. Lett. 79, 1 (1997); Q. A. Turchette et al., Phys. Rev. Lett. 81 3631 (1998). 
[9] D. Bouwmoostor, J-W. Pan, M. Daniell, H. Weinfurter, and A. Zeilinger, Phys. Rev. Lett. 82 1345 (1999); A. Rauschen- 
beutel et al. Science 288, 2024 (2000). 



6 



[10] C. S. Sackett, et al, Nature (London) 404, 256 (2000); J.-W. Pan, M. Daniell, S. Gasparoni, G. Weihs, and A. Zeilinger, 

Phys. Rev. Lett. 86, 4435 (2001). 
[11] W. H. Zurek, Phys. Today 44(10), 36 (1991). 

[12] D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore, and D. J. Jeomzem, Phys. Rev. A 46, R6797 (1992). 
[13] B. Julsgaard, et al, Nature 413, 400 (2001). 

[14] S. Mancini, V. Giovannetti, D. Vitali and P. Tombesi, Phys. Rev. Lett. 88, 120401 (2002). 
[15] S. B. Zheng, Australian J. of Phys. 53, 793 (2000). 
[16] S. Mancini and P. Tombesi, e-print quant-ph/0109138. 

[17] W. J. Munro, K. Nemoto, G. J. Milburn, S. L. Braunstein, e-print quant-ph/0109049. 

[18] A. Dorsel, J. D. McCullen, P. Meystre, E. Vignes, and H. Walther, Phys. Rev. Lett. 51, 1550 (1983); A. Gozzini, F. 

Maccarone, F. Mango, L Longo and S. Barbarino, J. Opt. See. Am. B 2, 1841 (1985). 
[19] L Tittonen, G. Breitenbach, T. Kalkbrenner, T. Miiller, R. Conradt, S. Schiller, E. Steinsland, N. Blanc, and N. F. de 

Rooij, Phys. Rev. A 59, 1038 (1999). 
[20] Y. Hadjar, P. F. Cohadon, C. G. Aminoff, M. Pinard, and A. Heidmann, Europhys. Lett. 47, 545 (1999); P. F. Cohadon, 

A. Heidmann and M. Pinard, Phys. Rev. Lett. 83, 3174 (1999). 
[21] A. Gillespie and F. Raab, Phys. Rev. D 52, 577 (1995). 
[22] M. Pinard, Y. Hadjar, and A. Heidmann, Eur. Phys. J. D 7, 107 (1999). 
[23] C. K. Law, Phys. Rev. A 51, 2537 (1995). 

[24] A. F. Pace, M. J. Collett, and D. F. Walls, Phys. Rev. A 47, 3173 (1993); C. Fabre, et al, Phys. Rev. A 49, 1337 (1994); 

K. Jacobs, P. Tombesi, M. J. Collett, and D. F. Walls, Phys. Rev. A 49, 1961 (1994); S. Mancini and P. Tombesi, Phys. 

Rev. A 49, 4055 (1994); G. J. Milburn, K. Jacobs, and D. F. Walls, Phys. Rev. A 50, 5256 (1994). 
[25] D. F. Walls and G. J. Milburn, Quantum Optics, (Springer, Berlin, 1994). 
[26] V. Giovannetti and D. Vitah, Phys. Rev. A 63, 023812 (2001). 

[27] Z. Y. Ou, S. F. Pereira, H. J. Kimble, and K. G. Peng, Phys. Rev. Lett. 68, 3663 (1992). 

[28] M. D. Reid and P. D. Drummond, Phys. Rev. Lett. 60, 2731 (1988); M. D. Reid, Phys. Rev. A 40, 913 (1989). 

[29] S. M. Tan, Phys. Rev. A 60, 2752 (1999). 

[30] M. D. Reid, e-print quant-ph/0112038. 

[31] L. M. Duan, et ai, Phys. Rev. Lett. 84, 2722 (2000); 

[32] R. Simon, Phys. Rev. Lett. 84, 2726 (2000). 

[33] V. Giovannnetti, S. Mancini, P. Tombesi, Europhys. Lett. 54, 559 (2001). 

[34] M. Pinard, P. F. Gohadon, T. Briant and A. Heidmann, Phys. Rev. A 63, 013808 (2000); T. Briant, P. F. Gohadon, M. 

Pinard, and A. Heidmann, e-print quant-ph/0207049. 
[35] A. Abramovici et al, Science 256, 325 (1992). 




FIG. 1. Schematic description of the system under study. L, being the equilibrium distance between the movable mirrors 
Ml, M2, is assumed to also be the distance between the fixed mirrors M3, Af4. The mirror M3 represents the input-output 
port of the cavity. 
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FIG. 3. Degree of entanglement E{uj) of Eq. (14) versus frequency and temperature T, in the case of a mechanical frequency 
mismatch ^2 — ^^i = 10 Hz. The plot has been cut at E{ui) — 1. The other parameter values are in the text. 
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